Second harmonic generation in metal nano-spheres: 
full-wave analytical solution with both 
local-surface and nonlocal-bulk nonlinear sources 
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We present a full-wave analytical solution for the problem of second-harmonic generation from 
spherical particles made of lossy centrosymmetric materials. Both the local-surface and nonlocal- 
bulk nonlinear sources are included in the generation process, under the undepleted-pump approxi- 
mation. The solution is derived in the framework of the Mie theory by expanding the pump field, the 
non-linear sources and the second-harmonic fields in series of spherical vector wave functions. We ap- 
ply the proposed solution to the second-harmonic generation properties of noble metal nano-spheres 
as function of the polarization, the pump wavelength and the particle size. This approach provides 
a rigorous methodology to understand second-order optical processes in metal nanoparticles, and to 
design novel nanoplasmonic devices in the nonlinear regime. 

PACS numbers: 42.65.Ky, 73.20.Mf, 78.35.-l-c, 78.67.Bf 



I. INTRODUCTION 

Nonlinear phenomena in metal nanostructures are 
gathering much attention due to their potential applica- 
tion as novel components for integrated optics More- 
over, second harmonic (SH) generation from metal nanos- 
tructures provides a powerful tool for probing physical 
and chemical properties of material surfacesi^i^ Noble 
metal nanoparticles support localized surface plasmons 
(LSPs). LSPs arc collective oscillations of the conduc- 
tion electrons, that strongly affect the optical response of 
the metal. When LSPs are resonantly excited, the local 
electromagnetic field is significantly enhanced in the par- 
ticle, enabling the generation of nonlinear optical effects, 
such as harmonic generation, at relatively low excitation 
powers. 

Second harmonic generation by nanoparticles origi- 
nates from two contributions: the particle bulk and 
the surface contribution, respectively. In noble metal 
nanoparticles, the local-bulk source is absent because of 
the material centro-symmetry, and only the nonlocal- 
bulk contribution needs to be consideredj^ The local- 
surface contribution to SH radiation is due to the sym- 
metry breaking at the interface with the embedding 
mediumi^i^ The magnitudes of the nonlocal-bulk and 
local-surface SH contributions depend on the shape of the 
nanoparticle and on the optical properties of the metal at 
the fundamental and second harmonic frequencies J^iii 

Both the bulk and surface SH sources contribute sig- 
nificantly to the SH scattering response of a spherical 
nanoparticle, as already demonstrated in approximated 
analytical models developed to describe the SH gener- 
ation from noble metal nanostructures with simple ge- 
ometries, such as spheres and cylinders. In 1999, Dadap 
et al. studied the SH radiation generated from the 



surface of a sphere with radius R much smaller than 
the wavelength of the incident light A {2ttR/X << 1, 
Rayleigh limit) It is shown that the leading-order con- 
tributions to SH radiation arise from the electric-dipole 
p(2i^) a^jjj ^i^g electric-quadrupole Q^^") moments, and 
the main selection rules for the SH scattering from a 
sphere are outlined. In Refs. [l^ and [l^, the Rayleigh 
limit is analyzed by taking into account both the bulk 
and the surfaces polarization sources, showing that the 
SH field is radiated by an effective electric dipole mo- 
ment pi/7^(f) = p'^"^ -h i fcoQ^^"^f/3 (the SH magnetic 
dipole emission is forbidden because of the axial sym- 
metry of the system). It is demonstrated that both the 
nonlocal-bulk and the local-surface SH sources contribute 
to the induced electric dipole moment p*^^"^', while only 
the local surface sources ^ntribute to the induced elec- 
tric quadrupolc moment Q'^"\ in the Rayleigh regime. 
The presence of distinct SH sources with their own radi- 
ation patterns causes the SH Rayleigh scattering process 
to differ significantly from the linear Rayleigh scattering. 
In particular, the theory predicts the absence of the SH 
signal in the forward direction and the (koR)^ scaling of 
the SH scattering cross-section. 

The SH Rayleigh scattering model is inaccurate if the 
particle size is comparable with the wavelength, owing to 
the contributions of SH multipolar orders higher than 2 
are not negligible. A full- wave analysis of the SH scat- 
tering from spherical particles of arbitrary size is devel- 
oped in Ref. [ij, but only the surface SH source was 
taken into account. Moreover, the enforced boundary 
conditions are incorrect resulting in zero SH radiation 
from the radial component of the source. Recently, a full 
wave theory of the SH radiation generated by a chain 
of parallel infinitely long cylinders, including both the 
bulk and surface nonlinear sources, has been developed 
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m Ref. i. A full-wave theory of the SH radiation gener- 
ated in three-dimensional structures consisting of metal- 
lic spheres made of centro-symmetric materials has been 
proposed in Ref. [isl but here again the treatment is 
limited to the surface source. 

In this paper we propose a full- wave analytical solution 
for the SH scattering from nano-sphercs of arbitrary size, 
which includes both the contributions of the nonlocal- 
bulk and the local-surface SH sources. As a case of con- 
siderable interest, we apply this solution to investigate 
the SH radiation generated by noble metal nano-spheres. 
The present work is organized as follows. In Section Hi Al 
the electromagnetic formulation of SH scattering by a 
metal nano-sphere is presented, and the analytical solu- 
tion for the fields is defined, at both the pump and SH 
frequencies. In Section IIIII the SH scattering from gold 
nano-spheres with increasing size is studied, as function 
of the wavelength and polarization of the pump field. In 
Section lTVl the conclusions are outlined. This manuscript 
is completed by six Appendices, where detailed formulas 
for the analytical calculation of all the quantities of in- 
terest are provided. 



II. PROBLEM STATEMENT AND SOLUTION 

A. Problem statement 

Let us consider the electromagnetic field at frequency 
2uj generated from a metal sphere of radius R, when illu- 
minated by a time-harmonic electromagnetic plane-wave 
at frequency oj incoming from infinity. We use a spher- 
ical coordinate system (O, r, 9, 0) with the origin O in 
the center of the sphere, as in Fig. [TJa); we denote 

with ^f, 0, ^ the unit vectors of the spherical coordi- 
nate system. The domain of the electromagnetic field 
is the entire space M'^, divided into the interior part of 
the metal domain , the embedding medium Oe and the 
metal surface S. The surface E is oriented in such a way 
that its normal n points outward, n = f |s . We use the 
convention a (r, = i?e {A*^^' (r) cxp (ifit)} for repre- 
senting a time harmonic electromagnetic field at angular 
frequency f2, where r = rf . The second harmonic gen- 
eration problem involves two electromagnetic fields oscil- 
lating at different frequencies: the fundamental-frequency 
electromagnetic field (E'^'^\ H^^^) at frequency w and the 
second harmonic electromagnetic field (E^^'^^, H^^")) at 

frequency 2a;. We denote with ^Eq"'', Hq"'^ the incident 
(pump) electromagnetic field: 

= (f ) (fco X so) e--r-r , 

where Eq is the electric field magnitude of the linearly- 
polarized pump beam, eg is its polarization direction. 




(a) (b) 



FIG. 1. (a) Scheme of the spherical particle and coordinate 
system {O, r, 9, (p). (b) Elementary closed curve Al across the 
selvedge region at the particle boundary. 

ko its propagation direction, kp"' = kofce('^), fce(w) = 
w^/ieEe and Ce = \/ fJ-e/^e- The parameters Eg and /ig are 
the permittivity and the permeability of the embedding 
medium. 

Since the intensities of the SH fields generated by noble 
metals are always orders of magnitude weaker than the 
intensities of the pump fields, the SH fields do not signifi- 
cantly couple back to the fundamental fields (undepleted- 
pump approximation). As a result, the electromagnetic 
scattering problems at the fundamental frequency and 
at the SH frequency are both linear. The linear elec- 
tromagnetic response of the metal is characterized, in 
the frequency domain, by the permittivity e^, which de- 
pends on the frequency, and by the permeability fii , that 
we assume independent of the frequency. The nonlinear 
optical response at SH frequency is phenomenologically 
described by a nonlinear polarization induced by the fun- 
damental electric field E*^"). This polarization acts as 
the source of the SH electromagnetic radiation, and has 
both a local surface contribution Pi^"^' defined on S and 
a nonlocal bulk contribution P^"'' defined in tli^^ 

Noble metals arc ccntrosymmctric materials, therefore 
the leading term of the bulk nonlinear polarization den- 

sity IS given by: 

Pf = eo xr^E(-) (VE(")) inli, , (2) 

„(2w) 

where Xt the quadrupolar contribution to the 

second-order nonlinear bulk susceptibility of the metal 
and Eo is the vacuum permittivity. Noble metals are also 
isotropic, therefore relation (2) becomes: 

Pr' =eo /3E("V-E(") + 

eo 7V (e(") • e(")) + eo 5' (e'") • v) E^") , 

(3) 

where 7 and 5' arc parameters depending on the ma- 
terial. Due to the homogeneity of the material, we also 
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have V • E^"^^ = 0, therefore the first term in ^ vanishes. 
Moreover, following Ref. [l^ we neglect the third term 
in (3). For these reasons we shall consider the following 
expression for the bulk SH source term: 



7V (: 



The surface nonlinear polarization contribution P 
given by: 



{2u 



(4) 



IS 



by gold spherical nanoparticles with R = ISOnm at 
A = 800nm, Bachelier et al. have found that an opti- 
mal choice for the phenomenological parameters a, b and 
d should be (a = 0.5 - iO.25, b = 0.1, d = 1)^ 

The order of magnitude of the bulk SH polarization is 
thus given by: 



5(2") 



Xb (w) Eo , 



(8) 



p(2t^) 



(2").eMe(") onS , 



(5) 



where Xs the second-order surface nonlinear suscep- 
tibility. Since the nanoparticle surface has an isotropic 
symmetry with a mirror plane perpendicular to it, the 

surface susceptibility tensor has only three non- 

vanishing and independent elements, X-LJ— L: X±|||| 
X||-L|l = where _L and || refer to the orthogonal and 

tangential components to the nanoparticle surface (Fig. 
(Ha)). Following Ref. [H we neglect the component x±\\ || j 
therefore: 



p1'"^ = eo [x±±±nhh+ 

Xiixii (tifiti+tants)] :E(")e('^) 



(6) 



In order to calculate the SH radiation generated by the 
metallic sphere, we have to evaluate: 

1. the electric field E'") in tli and on E, induced by 



the pump electromagnetic field ( E, 



^0 







2. the bulk and surface SH sources generated by E^'^^; 

3. the SH electromagnetic fields (E^^'^), H^^'^)) radi- 
ated by the bulk and surface nonlinear polarization 
fields. 

Both problems 1) and 3) are solved by expressing the un- 
known fields in terms of spherical vector wave functions 
(SVWFs), defined in Section [HB] In Section [iHl we will 
study the SH generation from gold nano-spheres, using 
as values for xJ-±J-i X||-L|| ^-i^cl 7 the expression in terms 
of the Rudnick-Stern (R-S) parameters {a,b,d), as done 
by Bachelier et al.)^^— 



X\\±\\ 
7 



Xb (w) 



p 



eo 



uj^ eno 
eno 



Xb (w) ■ 



eo 



eno 



(7a) 
(7b) 
(7c) 



where x& = (ej/eo — 1) is the bulk linear suscepti- 
bility of the metal, e is the electron charge (absolute 
value), no is the number density of the conduction elec- 
trons and Up is the plasma frequency. By choosing 
(a=l,6 = — l,ci=l), we obtain the Rudnick-Stern hy- 
drodynamic modelf^ By measuring the SH generated 



where Ap = cq (27r/a;p) and Et ~ (Ap/27r) eno/ eo- For no- 
ble metal nanoparticles Ap = 140 nm and « 10^*V/m. 
The order of magnitude of the surface SH source is in- 
stead given by Pi^"' = P^^'^^ (Ap/27r). 

It is worth noting that the proposed analytical solution 
of the SH generation from a sphere is independent of the 
expression of X-L-L-Li X||-L|| 7- 



B. Electromagnetic fields at pump frequency 

The electromagnetic field at the fundamental fre- 
quency is the solution of the Maxwell's equations: 



V X H, 



-iujSi (w) E 



n X H, 



H 



nx (E'^'-Efe) 



X Efc) = 




= n X H 

= n X eJ, 



H(w) 
, ^ sc 

fiCJEeEic' 







on E 



in 



(9a) 



(9b) 



(9c) 



where ( E,^"^ H^"^ 



Tt('^) 



denote the fields in fli and 
denote the scattered fields in Qe, namely 



e1") = Ei") - E^") and Hfc^ = ui"^ - U^"^ 



Q txiiia j-j-sc — j-j-r — j-iQ . Equa- 
tions (O have to be solved with the radiation condition 
at infinity for the scattered fields. Due to the symmetry 
of the problem, the general solution of the source-free 
Maxwell equations are expressed in each homogeneous 
region through the SVWFs M^j 



'^^^ and nL^1)^>21 



Mlil(kr,0,. 



(kr) X„ 



N(;f)(A:r,0,0)--VxM(„^), 



k 



(10a) 



(10b) 



where z. 



(J) 



(kr) is one of the four kinds of the 



spherical Bessel functions, namely, Bessel function of the 
first kind jn = jn (kr), or Bessel function of the sec- 
ond kind y„ = i/n (kr), or Bessel function of the third 
kind (spherical Hankel functions of the first and sec- 
ond kind), h^P = hn '' (kr) and hn^ — hn^ {kr), which 
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we denote, respectively, with the apices J = 1,2,3,4; 
Xmri = (^'i 4>) IS & Spherical vector harmonic (Ap- 
pendix [A|. Here, both the spherical vector functions and 
the vector spherical harmonics are indexed by the order 
m and the degree n. 



The incident plane- wave ^Eq'^'', Hq"'^ is decomposed 

in Eq. (|lla[) through the regular SVWFs, non-singular 
in the center of the sphere (J — 1), where Eq is the 



amplitude of Eq, and the coefficients Ipmnj'ZmnI are 
given in Appendix B for a linearly polarized state along 
the x— axis. Also the unknown fields m 
r^i (0 < r < R), are decomposed through the regu- 
lar SVWFs in Eq. ([TTb| . where Ci (w) = (w), 
fcj (w) = u:^£i (bj) Ui- The unknown scattered fields 



Mi 

I^eIc'', hIc'' j in f2e ( for i? < r) are instead decomposed 

in Eq. (|llc[) through the radiative SVWFs, satisfying the 
radiation condition at infinity (J = 3). 



oo n 



E, 



n—1 m— — n 
j-^ oo n 



n—1 m— — n 



(lla) 



cxD n 



Ei"' (r. 



Ti— 1 m— — n 

oo 71 

r ( \ {dinlMl^l [h (w) r, e, 0] + ci-,Xl [h (co) r, 0, 0] } , 



En 



(lib) 



. m— — n 



Efe' (r,0, 
Hi^f) (r,0,0) 



OO n 



^oE E {^^^2M(^Ufce(a.)r,0,<^]-f a(;rlN(^U^e(u;)r,0,,/)]} 

n—1 m— — n 

oo n 

-T^E E (a^lMW [fce(w)r,0,0]+6(;r2N« [fce(c.)r,0,0]}. 



(11c) 



n—1 m— — n 



E^e^mcM satisfy Eq. 



The decomposition (HTc]) of (^Eic\ H^' 
(|9cl) . and the decomposition (fTTb]) of (^e|'^\ h|'^^^ sat- 
isfy Eq. (Pa|) . The unknown coefficients |am,l,6m2,| 

and |cm2,(imri| are determined by requiring that the 
decompositions (|llblllcp also satisfy the boundary con- 
ditions (|9bp . The analytical expressions of |amTi,&rml| 

and |cTO2,rfTrm| as functions of {Pmn^Qmn} are given in 
Appendix [Cl 



C. Electromagnetic fields at second harmonic 

The SH electromagnetic field satisfies the Maxwell's 
equations: 



V X E 

V X H 



(2a;) ^ 
(2w) 



(2") 



2iuje^ {2uj) Ef -t- J 



J2u 
'b 



in fii, 
(12a) 



n X ( H.^ 
n X ( E, 



H 



(2u 



E, 



{2u 



_ _,(2c^) 
Jelet 

,(2m) 
■ jmag 



V X E, 



V X H, 



{2uj) 



(2") 



on S 



where 



:(2") 
ielet 



— — 2iwn X ( n X P 



>(2": 



(12b) 



(12c) 



(13a) 



(13b) 



(13c) 



the operator denotes the surface gradient, 

(Ef"\Hf"') denote the SH fields in i7„ (Ei'^^Hi""^) 
denote the SH fields in 51^ and e' is the selvedge region 
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(see Appendix |D]) pcrmittivityji^ which we assumed 
equal to Eq. The sources of the SH radiation, therefore, 
are of three types. The volume current density field 
J^^""* given by Eq. (|13ap . takes into account the con- 
tribution of the SH bulk nonlinear polarization. The 
surface electric ciurrent density given by Eq. (jl3bp . 
takes into account the contribution of the SH tangent 
surface nonlinear polarization. The surface magnetic 

current density jmag given by ()13c|) . takes into account 
the contribution of the SH normal surface nonlinear 
polarization (see Appendix D). The systems of Eq. (|12ap 
have to be solved with the radiation condition at infinity 

for (Ei^"\Hi^'^)). 

The SH field equations are formally the same of the 
fundamental field equations except for the bulk source 
term in the Maxwell- Ampere equation and the 

substitutions w — >■ 2a;. Consequently, the problem is 
reduced to that already solved for the fundamental fields, 

I 



by expressing the electric field inside the nanoparticle 
as: 



where 



E 



E 



{2u:} 
part 



(2m) 



E 



{2u 



E, 



(2cL;in n„ 
part 



(14) 



Ei (2a;' 



-7V (e(") • E(^) j . (15) 



The fields (E(t:!,Hf")) and (Ei''^\ H^'")) are solu- 

tions of the homogeneous Maxwell equations at frequency 
2a; and have to satisfy the boundary equations: 



f X H 



(2a;) 



H 



(2a;) 



„:(2") 
Jelet 



— Jmag 1 A XL; 



(2a;) 
part 



on S 



Therefore we impose: 



(16) 



Hf-) (r,0,, 



oo n 



E^t! ir, e, 0) = ~Ei'-^ E [fc^ (2..) r, 9, cj,] + ^'^^^nW [fe^ (2a.) r, 0, cj,] } 



n—1 m——n 
(2cj) oo n 



iQ (2a;) ^ ^ 

^ ^ n—1 m——n 



for r < and (17a) 



^ ^ {dl^^'>Ml^l [h r, 0, 0] + c(2^)N(£ [h (2c.) r, 9, cf\] 



OO 71 

E(2") (r, 0, ,/)) = ii;^") ^ 5] {?>^^^)M(^), [fc, (2a;) r, 0, 0] + a(^:^)N(^), [K (2a;) r, 0, 0] } 



H(2") (r,0. 



n—1 m— — n 



(2a;) 



OO n 

^ 5^ {a(^-)MW [fc42a.)r,0,0]+C:^)NW [fc42a.)r,0,^]} 



for r > R, 



(17b) 



n—1 m— — n 



r 



where 



E. 



(2a;) 



Xb (lo) Eo (18) 



is a characteristic electric field expressing the order 
of magnitude of the SH electric field in the nanopar- 
ticle. Both the expressions of ^e|j^^"j|, Hp"-*^ and 

(Ei^"\Hi^"^) satisfy the homo geneous Maxwell equa- 
tions at frequency 2aj. The unknown coefficients 
I 



^a^n\b^n^^ and |c^^5^\ (imn''| are evaluated by impos- 
ing that the functions ([T7)) satisfies the boundary equa- 
tions (|T6|) . The contribution of the SH bulk source, there- 
fore, is given by the term —fx Ep'^") |s in the system 
of Eq. (|16p . playing the same role as the surface mag- 
netic current jmag ■ In order to analytically evaluate the 
I and |c^^!;^\ the fields 

- f X E^^„"J (r = R, 9, (/)) are ex- 
panded in terms of SVWFs: 



and i man 



coefficients \ amm\ fe*'^"' 

.(2a;; 
Jelet 



,(2a;) 
Jelet 



;(2a;) 



E, 



(2a;) 



OO n 

E E p 

n—1 m— — n 



(2a;)„ 
mn -^rnn 



+ " mri ^ ^jriTi 



f xE^tl (r = i?,( 



_^(2.) ^ 

n—1 m— — n 



n 

E 



(19) 



(2a;)- r- Y 
mn r X [r X ^mn 



{OA)] 



where the coefficients |t;'m^' , u'mSl^'' | and are evaluated in Appendix E. The formulas of the coefficients 



(2a;) 



and Vmn^ obtained in this way arc given in Appendix iFl 
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FIG. 2. SH scattering geometry. The pump electromagnetic 
field propagates along the positive direction of the 2— axis 
and is linearly polarized along so- The SH scattered field 
is observed along the direction K = x (The scattering plane 
is xOz). The y— and z— components of the SH field are 
considered, respectively parallel (||) and orthogonal (_L) with 
respect to the scattering plane. 

III. DISCUSSION: GOLD NANO-SPHERE 



the origin of the coordinate system. — ^'dn — depends 
on the collection direction K of the scattered SH light. 
cica^ has the physical dimensions of an area, and is pro- 
portional to the SH generation efficiency. In order to 
analyze the SH radiation polarization state, the analyzer 
can be polarized either parallel (|j) or perpendicular (_L) 
to the SH scattering plane, defined by the propagation 
direction feo of the pump wave, and the collection direc- 
tion (K). We denote with dPl^^'^''/dn and dP^^"'>/dn the 
radiated powers per unit solid angle associated to the || 
and _L components. The analysis of the polarization state 
of SH radiation collected at right angle from the pump 
beam, i.e. K = x, is very important because it allows 
to discriminate the radiation generated by even and odd 
SH multipole sources. Indeed, only the SH N„i„ multi- 
poles with odd n contribute to the || component and only 
the SH Nill multipoles with even n contribute to the _L 
component. This behavior is reversed for the SH Mmn 
multipoles. 



In the present section, by using the analytical solu- 
tion derived in the previously, we analyze the SH gen- 
eration from an isolated gold nano-sphere in vacuum, 
as the radius, the pump wavelength and polarization 
vary. Specifically, we study the SH radiation generated 
at pump wavelengths of A = 780 nm (Ti:sapphire laser) 
and A = 520 nm (gold plasmon resonance). Particular 
care has been devoted to the comparison with the exist- 
ing theories in the Rayleigh regime. In order to model the 
bulk linear susceptibility of gold, we interpolated John- 
son and Christy's experimental datai^ We used two sets 
of values for the R-S parameters (see Eq. ^) of the 
SH sources: the set (a = 1, 6 = — 1, d = I) proposed in 
Ref. [H and the set (a = 0.5 - i 0.25, b = 0.1, d = 1) 
proposed in Ref. [13. 

The pump electromagnetic field is a plane- wave prop- 
agating along the positive direction of the z— axis, and 
linearly polarized in the xy plane, with a polarization di- 
rection fo- We indicate with a the angle between the unit 
vector £o smd x— axis (the reference versus is counter- 
clockwise, seen from the semispace z > 0), as shown in 
Fig.H 

We investigate the SH power per unit solid angle 

— — -, radiated in the farfield along the direction K 
and collected by an analyzer with polarization state e*: 



dPi?"^ ( K 



dn 



lim 

r— ^oo 



2Ce 



■g(2aj; 



(20) 



and the SH scattering cross-section cica 

2 

Cil-J = fim 



E. 



E, 



(21) 



where is a spherical surface with radius p, centered at 



A. SH source currents: Rayleigh and Mie regimes 

Here we analyze the SH radiation generated from the 
single nonlinear sources, acting as if they were radiat- 
ing independently. Fig. [3] shows the magnitude of each 
SH source current density, namely jf"-* , jfi'^^, jma], 
normalized to its own maximum, and computed for the 
two pump wavelengths A = 520 nm and A = 780 nm, 
corresponding to resonance and off-resonance conditions, 
respectively. Two nanoparticle radii have been consid- 
ered, namely a particle with small size (i? ~ 10 nm) and 
a particle comparable in size to the pump wavelength 
(R ~ 150 nm). The pump field is linearly polarized along 
the X— axis (a = 0). In particular, the first column (pan- 
els a-d) shows the magnitude of the bulk electric current 
density J^^'^'' in the xOz plane, while the second and 
third columns (panels e-h,i-l), show the magnitude of the 

surface electric jg;^"* and magnetic ]mag current densi- 
ties, respectively. For small particles, is significant 
across the entire particle volume (panel a, b). In par- 
ticular, while for A = 520 nm (panel a) J^^"^-* decreases 
along the direction of forward scattering, for A = 780nm 
it is almost uniform. For particles with larger size (pan- 
els c,d), the skin effect appears, i.e. the current J^"'' is 
strongly confined near the particle surface. The intensity 
distribution of both jg^^j' and jmag (c,f,i,j) is symmetric 
around the polarization direction of the pump field for 
small particles. As the radius increases, this holds no 
longer true due to the onset of higher order multipoles 
(panels g,h,k,l). It is worth noting that the surface elec- 
tric current density jg;"^'' vanishes on a circle lying in the 
yOz plane, for any particle size and pump wavelength, 
as shown in panels (e-h) , due to the rotational symmetry 
of the particle. Similarly, the surface magnetic current 
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FIG. 3. SH source current distribution excited by a pump field 
linearly polarized along x. The nano-sphere size is i? = 10 nm 
(first and second row)and R = 150 nm (third and fourth row), 
the pump wavelength is A = 520 nm (first and third rows) and 
A = 780 nm (second and fourth rows). Panels (a,b,c,d) are 
relative to the bulk current density cut in the xOz plane, pan- 
els (e,f,g,h) are relative to the surface electric current density 
and panels (i,j,k,l) are relative to the surface magnetic current 
density. Each panel shows the current magnitude normalized 
to its own maximum. 




FIG. 4. SH power per unit solid angle collected at right an- 
gle from the pump beam (Fig. ^ as function of the pump 
polarization angle a, for a nano-sphere with R — 10 nm. The 
pump wavelength is 780 nm. The blue line corresponds to 
the II component and the red line to the _L component with 
respect to the scattering plane. All the graphs are normalized 
to the maximum of the most intense component. Panel (a) is 
relative to the bulk current, panel (b) to the surface electric 
current, panel (c) to the the surface magnetic current, as if 
they acted separately. 



density jmag displays a circle with a constant magnitude 
for every particle size, as shown in panels (i-1). 

In Fig. mthe SH power per unit solid angle collected at 
right angle from the pump direction is shown for a nano- 
sphere with R — 10 nm. Panels (a), (b) and (c) are rel- 
evant to the SH radiation generated by the bulk, surface 
electric and surface magnetic SH source currents, respec- 



tively, as if they acted separately. The blue and red lines 
correspond respectively to dPy^"-'/'^^ ^iid dP^"^ /dQ , 
and for each panel both the components are normalized 
to the maximum of the most intense. 

These results agree with those obtained analyzing the 
SH radiation from a metal nano-sphere in the Rayleigh 
limiti^ii^ In this regime, the SH radiation coincides with 
the electromagnetic field radiated by a fictitious elec- 
tric dipole with effective moment p^^^j}{r) ^ p(2'^) -|- 

i fcoQ'^^"'f/3, where p^^") is the induced SH electric 

dipole moment (71 = 1) and Q^^"^) is the induced SH 
electric quadrupole moment (n = 2). Depending on the 
component of the SH intensity, two different shapes of 
the polarization diagrams arise. For each SH source, 
rfPip^'/dri is due to a SH dipolar electric mode aligned 
along the propagation direction of the pump^^ii^ there- 
fore its value is independent of the polarization angle of 
the pump. On the other hand, the four lobe pattern ob- 
served for dPf"'^ /dO, is due to a SH quadrupolar mode. 
Furthermore, dP^^^ / dfl is negligible for both and 
je?et^ while it is comparable with dP^^^'^'^ / dfl for jmag. 
For larger particles, higher order SH multipolcs arise due 
to larger retardation effects, significantly modifying the 
SH radiation characteristics, as we shall see in the next 
section. 



B. SH scattering cross-section 

In this Section, we study the SH scattering cross- 
section cica ^ for a gold nano-sphere, using as R-S param- 
eters the values {a = l,b = —l,d = 1). All the results are 
relative to the case of a pump plane-wave with electric 
field of unitary magnitude, i.e. \Eo\ = 1 Vm~^ . Figure 
[5] shows cica^ as function of the pump wavelength (black 
lines), for four particle sizes. For all the investigated 

sizes, Csca"^ shows a maximum at A ~ 520 nm, when the 
pump wavelength matches the plasmonic resonance of 
the particle. Another local maximum is also observed at 
A sa 1040 nm. At this wavelength the SH fields resonate 
in the gold nano-sphere. The relative intensity of cica"^ 
at A « 1040 nm increases as the particle size increases. A 
third local maximum can be observed at A « 700 nm, for 
certain particle sizes {e.g. R = 100 nm and R ~ 200 nm). 
Similar trends, not shown here, have been found with the 
set of values (a = 0.5 - zO.25, 6 = 0.1, d = 1) for the R-S 
parameters. 

In order to unveil the origin of the SH radiation in the 
Rudnick-Stcrn model, the contributions of each multipole 
are shown, up to the 6th order. For R = 10 rtm (panel 

a), Csca'' is mostly due to the SH dipolar source, and 
only for short wavelengths the SH quadrupolar source 
begins to be significant. As we increase the radius R to 
100 nm (panel b), we can identify three different regimes: 
for short wavelengths A < 550 nm, G^a'^ is dominated 
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FIG. 5. SH scattering cross-section (black line) as function 
of the pump wavelength for nano-spheres with R — 10 nm 
(a), J? = 100 nm (b), R = 150 nm (c), R = 200 nm (d), with 
{a = l,b = —1, d = 1). The contribution of each multipolar 
order to the total radiation cross-section is shown up to the 
6th: n = 1 (green), n = 2 (blue), n — 3 (violet), n = 4 (cyan), 
n = 5 (red), n = 6 (yellow). 
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FIG. 6. SH scattering cross-section as function of the nano- 
sphere radius R at pump wavelengths A — 520 nm (red), 
A — 780 nm (green), A = 1040 nm (blue), obtained with 
(a= l,fe = -l,d= 1). 



by the octupolar source, the quadrupolar one prevails 
in the range 550 nm < A < 950 nm, while the dipolar 
source is the most intense for large wavelengths. For a 
particle with R = 150 nm (panel c), the dipolar source 
is negligible regardless of the pump wavelength, and the 
main contributions to cica"^ arise from the multipoles 
with 71 = 2,3,4. Similarly, the main contributions for a 
particle with R = 200 nm arise from n = 2,3,4 for large 
wavelength, and from the multipoles n = 5, 6 for short 
wavelengths. 

Figure ini shows the SH scattering cross-section as func- 
tion of the nanoparticlc radius, for three values of the 





FIG. 7. SH power per unit solid angle collected at right 
angle from the pump beam (Fig. [2]) as function of the 
pump polarization angle a, for nano-sphere with R = 10 nm 
(a,e), R = 100 nm (b,f), R = 150 nm (c,g), R = 200 nm 
(d,h), obtained by using (a-d) (a = 1,6 = —l,d = 1) and (e- 
h) (a = 0.5 — iO.25, b — 0.1, d — 1). The pump wavelength is 
A = 780nm. The blue line corresponds to the || component 
and the red line to the 1. component with respect to the scat- 
tering plane. All the graphs are normalized to the maximum 
of the most intense component. 



pump wavelength (A = 520 nm, 780 nm, 1040 nm). The 
SH scattering cross-section increases more than 4 or- 
ders of magnitude when the particle size increases up 
to i? = 100 nm. For larger radii, the SH scattering cross- 
section saturates and a small modulation takes place. 
For small particle size, the highest cross-section is shown 
when the particle is in plasmonic resonance (i.e., red 

curve). For larger particle size, the magnitude of cica'* 
is comparable for all the investigated pump wavelengths. 
Also in this case, similar trends have been found using the 
set of R-S parameters [a =- 0.5 — zO.25, h = 0.1, d = 1). 



C. SH power dependence on the pump polarization 

The SH power radiated at right angle from the prop- 
agation direction of the pump allows for the recognition 
of even- and odd-order multipolar contributions to the 
SH generation process, as already pointed out. Fig- 
ure [7| shows (iP||^'^''/(iri and dPf'^^/dfl as function of 
the polarization angle a of the pump, for two different 
choices of R-S parameters. The first row (a-d) is rel- 
ative to the R-S parameter set {a = l,b = —1, d = 1), 
while the second row (c-h) is relative to the set 
(a = 0.5 — zO.25, b = 0.1, d = 1). Four different values of 
particle size are presented: R = 10 nm (a,e), R = 100 nm 
(b,f), R = 150 nm (c,g), R = 200 nm (d,h). Both of 



dPl^'^^/dn and dP]_"'"'/dn are normalized to the max- 
imum of the most intense component. For small (a,e) 
and very large (d,h) radii, dPp^'/drj prevails over the 



dP'f^' /dfl for both sets of the R-S parameters. For lu- 



ll 



termediate sizes (b,c,f,g), dPj^^'^^/dn and dPf'^'/dn are 
comparable, and their relative intensities strongly depend 
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FIG. 8. SH power per unit solid angle of the || component 
collected at right angle from the pump beam (Fig. [2]), as 
function of the pump polarization angle a, for nano-spheres 
of size R = 80 nm (green), R — 120 nm (red), R = 140 nm 
(black), R = 150nm(blue) and (a = 1,6 = = 1). The 

pump wavelength is A = 780nm, all the curves are normalized 
to their own maximum. 



on the particular choice of the R-S parameters. For small 
particles, dP|p"''/(i51 is independent of the polarization 
angle a (panels a,e), while up to 4 lobes can appear for 
larger particles, as it will be shown more in detail in Fig. 
m The graphs of the component dPf'^'' /dfl feature four 
lobes oriented along the bisectors of the 4 quadrants, in 
each of the investigated case. 

The appearance of an octupolar SH source significantly 
modifies the shape of c?P|p"''/'^^- the Rayleigh limit 
the II component generated by the SH dipolar source 
fully prevails over the _L component generated by the 
SH quadrupolar source. As the radius increases, the in- 
tensity of the SH octupolar source increases due to the 
retardation effects, as pointed out in the previous section. 
Due to the interference in the far field of the SH dipolar 
and octupolar fields, the dP^'^'^'' / dfl reduces significantly 

until it becomes smaller than dP^^'^^ / dVL (b,c). As the ra- 
dius further increases, the SH octupolar source prevails 
over the SH dipolar source, and the shape of dPlj^'^^ /d^ 
gets close to that of an octupole (d,e). 

The details of the transition from the dipole to the oc- 
tupole pattern as the radius increases are shown in Fig. 
m First, the circular shape of dP'^^'^^ /dO. is shrunk along 
the directions a = 0°, 180°, until the amplitude reaches 
a null (green and red curves). Then, two lobes arise along 
these directions (black curve), forming a four- lobe pat- 
tern. As the radius further increases, the intensities of 
the two lobes along the directions at a = 0°, 180° pre- 
vail over the intensities of the lobes along the directions 
at a = 90°, 270° (blue curve). In conclusion, even some 
discrepancies arise from the comparison between the re- 
sults obtained by using the two sets of R-S parameters, 
the main features of the SH radiation field are common. 




FIG. 9. SH radiation diagrams for nano-spheres of size 
R = 10 nm (panels a,b), R = 150 nm (panel c, d), obtained 
by using (a = 1, 6 = — 1, d = 1). Panels (a, c) are relative 
to the pump wavelength A = 520 nm, and panels (b,d) to 
A = 780 nm. All the intensities are normalized to their own 
maximum. 

D. SH radiation diagrams 

Figure [9] shows the angular distribution of the SH radi- 
ation generated by gold nano-spheres. obtained using the 
R-S parameters (a = 1, b = —1, d = 1). The first row 
(a,b) is relative to a small sphere {R = 10 nm) and the 
second row (c,d) to a large sphere {R = 150 nm). In the 
first column (a,c) the pump wavelength is A = 520 nm, 
while in the second column (b,d) A = 780 nm. 

For particles with R = 10 nm, the dipolar and 
quadrupolar SH sources dominates the response, in 
agreement with the Rayleigh limit. In particular, we no- 
tice that the quadrupolar SH source is more important 
at A = 520 nm, while at A = 780 nm the dipolar mode 
fully dominates the response. As the particle radius in- 
creases, higher order modes come into play, resulting in 
an higher number of secondary lobes. For particles with 
R ~ 150 71m, the octupolar mode dominates the angular 
distribution of the SH radiation at A = 780 nm. More- 
over, for each particle size and pump wavelength, the 
lobes display a preferential alignment along the polariza- 
tion direction of the pump field. 

Figure [TO] shows the two cuts of the 3D radiation dia- 
gram along the xOz, and yOz plane. It is worth noting 
that the SH power vanishes in both the forward- and 
backward-scattering directions, regardless of R. This se- 
lection rule is a direct consequence of the rotational sym- 
metry of the sphere around the propagation direction 
of the pump. Two lobes, directed orthogonally to the 
pump propagation direction, characterize the radiation 
diagrams of small particles. As the particle size increases, 
the lobes with higher power tend to come closer to the 
forward direction, while the lobes in the backward direc- 
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FIG. 10. SH radiation diagrams as function of 9 angle, at 
tp = Q° (blue) and ip — 90° (red), for nano-spheres of size 
R = lOnm (a), R = 100 nm (b), R = 150 nm (c), and R = 
200 nm (d) and (a = 1, fc = —1, d — 1). The pump wavelength 
is A = 780 nm, for each panel the curves are normalized to 
the maximum of the most intense one. 



tion decrease in amplitude. The same trend is observed 
if the particle size is fixed and the pump wavelength de- 
creases. This behavior has been already observed exper- 
imentally for silver nanoparticles in Ref. [23l 



where 



' mn 



(cose) =_p-(cos^ 

(cose) ^±p^[cose) 



and P™ = P™ (u) is the associated Legendrc function of 
the first kind and of degree n and m. 



Appendix B: Calculation of |pmn,gm2| 



The expansion coefficients in Eq. (jllap . for a linearly 
polarized plane- wave propagating along the z— axis with 
the electric field parallel to the a;— axis (Fig. [1^), are: 



Pin = lin = -P-1 



IV. CONCLUSIONS 

We have developed a full-wave analytical solution for 
the second-harmonic generation from a spherical particle 
of arbitrary size, made of lossy centrosymmetric material. 
This approach extends the existing theories, enabling a 
rigorous treatment of all the sources of SH radiation, lo- 
cated both on the surface and in the bulk of the particle. 
The solution of the problem is derived in the framework 
of the Mie theory by expanding the pump field, the non- 
linear polarization sources and the second harmonic fields 
in series of spherical vector wave functions, and enforcing 
the boundary conditions at the sphere surface. 

Using the developed method, we investigated the 
spatial distribution of bulk and superficial non-linear 
sources. We also studied the SH radiated power as func- 
tion of the pump polarization angle, the SH radiation 
cross section spectrum and the radiation diagrams. 

The application of the proposed method in combina- 
tion with experimental observation, can improve the gen- 
eral understanding of nonlinear processes in metals, and 
can lead to an accurate evaluation of weights for the 
different SH sources. The theory of SH scattering can 
be easily extended to the multiparticle case. This ap- 
proach can also guide the design of novel nanoplasmonic 
devices with enhanced SHG emission for a wide range of 
applications;^ including sensors for probing physical and 
chemical properties of material surfaces. 



Appendix A: Vector spherical harmonics 

The vector spherical harmonics Xm„ arei ^^i^^ 
1 2n+l{n-m)\ 



^Jn (n + 1) y 47r [n + m)\ 

ilTrnn (cOS 9)9- Tmn (cOS 9) e'""^ , 



Appendix C: Calculation of |am2,6™|, |cm2,d™| 



The coefficients 
pressed as 



{UJ) 

0-rnn 

IP) 

Pmn 



■7 Li_ 



Pmn 



where Xe 



1^) 



and V'n = "ipn (p) , 



in = in (p) are the Riccati-Bessel functions defined as 

V'n(p) = P in{p), in = P hn \p). C deuotcs the first 
derivative of ^ = (^{p) with respect to p. 



Appendix D: Selvedge region 

The selvedge region (Fig. [T})) is a layer of infinites- 
imal depth S at the interface metal-vacuum. In this 
region there is a volumetric current density J^^' = 
i2aj ^Pi^"' • n/(5j n, which is exactly compensated by 
the normal component of the displacement current den- 
sity, J^^^ 4- i2a;D'^'^'' = 0, otherwise there would 
be an unbounded magnetic field. Therefore, in the 

selvedge region, D^,^'^-' = — ^P^^"-* • n. From the 
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Faraday-Neumann's law, applied to the elementary curve 
shown in Figure lb, we have ^Ef"^' - E^^"^) |s • 
A^ii ^ U2 — Ml, where = /a;*"' E'^") • dl ~ 



(Pr' ■ ft) 



Q(») 



and a = 1,2. By combining these rela- 



tions we obtain the equations ft x ( E, 



ft X V. 



Appendix E: Calculation of |am5^\&m"'|, |c$,if;^', d^"' | 

The coefScients |amn\6mn''| and |cmn\'imn^| are 
expressed as: 

{2uj) ^ /(2c^) ,/(2l^) U2u) ^ ,„(2c^) 

inn rtin ' mn i rnn ?nn ~ ^ mn i 

„(2(^) ^ ^t(2u^) j_ ^//(2c^) ,(2t^) ^ ,/(2cu) ,„(2cu) 



Tnn ' mn 



where with one apex we denote the contribution due to 
the tangential surface SH sources and with two apices we 
denote the contributions of both the normal surface SH 
sources and the bulk SH sources. For the contribution of 
the tangential surface SH sources we have: 



„/(2c^) 

mn 



„(2")n 



,(2")^ 



J2c^)x 



!,/(2") 



(2t^). 



,,/(2<^) 

^ mn 




)t,nyXe ) 


„/(2") 
mn 








,,/(2") 
^ mn 




l'^(2'^)\ Ce ;e ^„(2'^) 

'^■^'^ Ci(2tj)^"^^ 




mn 









,/(2'^) 



ci^V'n(xf-))^„(x^)) - uxr>)MxY'^') 



where the coefficients are given in Ap- 

pendixE xi'") ^ kf'^^R, = A:f "^i?, ^„ = V« (p) , 

= (p) are the Riccati-Bessel functions. 
For the contribution of both the normal surface SH 

sources and the bulk SH sources we have: 



„'/(2") 

mn 

,//(2") 

mn 

,„(2(^) _ 
^ mn 

„"(2") , 
mn 

J//(2'^) 

^ mn 

,//(2") 



Ci(2i.i;)'^'= Vn\Xi ) 



{2u 



m„ ,w ,„„ I are given m Ap- 



where the coefficients 
pendix iFl 



Appendix F: Calculation of , w'm!;^' | , |?t"mn\ f"™' | 

The coefficients lu'^^', f'^^' | for the surface tangential source can be expressed as: 



,{2lu) _ „C 



ctt ni OO 712 



„/(2'^) 
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/- / J / J / J / J mini m2n2 nimi n2m2nm ' mini m2n2 nimin2n 

SO — 



m mi— — ni n2 m2 — — n2 
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where Cj^l^^j^f^j^ are the Clebsch-Gordan coefficients [Chapter 8 in Ref. [21|, and the quantities in braces are 
Wigner 6j and 9j symbols [Chapters 9 and 10 in Ref. |2l|. 

The coefficients |M"mri\ ''^"mn'' | both the bulk and the surface normal polarization source can be expressed as: 



i,(2La) . /—, — — p-g,„„(fc^(w)i?) eo . /—, — — rT./'mn(fcj(w)i?) 

U 'mn = «Vn n + 1) , I + — 7T-^«V'^ + 1) TTT^ 

ko[uj)R ei[2uj) ko[uj)n 
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